Norges Teknisk-naturvitenskapelige Universitet Efficient Time-symmetric Simulation of Torqued Rigid Bodies Using Jacobi Elliptic Functions
نویسندگان
چکیده
If the three moments of inertia are different from each other, the solution to the free rigid body (FRB) equations of motion is given in terms of Jacobi elliptic functions. Using the Arithmetic-Geometric mean algorithm, [1], these functions can be calculated efficiently and accurately. The overall approach yields a faster and more accurate numerical solution to the FRB equations compared to standard numerical ODE and symplectic solvers. This approach performs well also for mass asymmetric rigid bodies. In this paper we consider the case of rigid bodies subject to external forces. We consider a strategy similar to the symplectic splitting method proposed in [16]. The method here proposed is time-symmetric. We decompose the vector field of our problem in a FRB problem and another completely integrable vector field. In our experiments we observe that the overall numerical solution benefits greatly from the very accurate solution of the FRB problem. We apply the method to the simulation of artificial satellite attitude dynamics.
منابع مشابه
Norges Teknisk-naturvitenskapelige Universitet Fakultet for Informasjonsteknologi, Matematikk Og Elektroteknikk Hovedoppgave
.........................................................................................................................I PREFACE........................................................................................................................... II TABLE OF CONTENTS................................................................................................... III
متن کاملNORGES TEKNISK-NATURVITENSKAPELIGE UNIVERSITET A method of parametric solution of convolution equations
A variant of the method of moments is developed for parametric solution of convolution equations of the first kind. Two models — the gamma model and the shifted gamma model — are studied in details.
متن کاملNORGES TEKNISK-NATURVITENSKAPELIGE UNIVERSITET Parameter Estimation in High Dimensional Gaussian Distributions
In order to compute the log-likelihood for high dimensional Gaussian models, it is necessary to compute the determinant of the large, sparse, symmetric positive definite precision matrix. Traditional methods for evaluating the log-likelihood, which are typically based on Choleksy factorisations, are not feasible for very large models due to the massive memory requirements. We present a novel ap...
متن کاملNORGES TEKNISK-NATURVITENSKAPELIGE UNIVERSITET Implementing Approximate Bayesian Inference using Integrated Nested Laplace Approximation: a manual for the inla program
This manual describes the inla program, a new instrument which allows the user to easily perform approximate Bayesian inference using integrated nested Laplace approximation (INLA). We describe the set of models which can be solved by the inla program and provide a series of worked out examples illustrating its usage in details. Appendix A contains a reference manual for the inla program. This ...
متن کاملNorges Teknisk-naturvitenskapelige Universitet Control Variates for the Metropolis-hastings Algorithm Control Variates for the Metropolis-hastings Algorithm
We propose new control variates for variance reduction in the Metropolis–Hastings algorithm. We use variates that are functions of both the current state of the Markov chain and the proposed new state. This enable us to specify control variates which have known mean values for general target and proposal distributions. We develop the ideas for both the standard Metropolis–Hastings algorithm and...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
دوره شماره
صفحات -
تاریخ انتشار 2006